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УДК 539.3+624.131 


Осесимметричная контактная задача консолидации 
для непрерывно-неоднородного по глубине полупространства! 


Л. Н. Евич 
(Донской государственный технический университет) 


Рассматривается постановка осесимметричной задачи консолидации для пористого, неоднородного по глу- 
бине полупространства и построение фундаментальног о решения для определения полей перемещений, 
деформаций, напряжений и порового давления при заданных граничных условиях. Отдельно рассмотрено 
решение уравнений описывающих напряжённое состояние порист. ОЙ среды под воздействием касательного 
усилия и под воздействием нормальной И радиальной нагрузках. При решении используется интегральное 
преобразование Ханкеля, которое позволяет свести задачу к решению системы обыкновенных дифференци- 
альных уравнений второго порядка. С помощью метода моделирующих Функций получены представления 
Для напряжений, смещений, порового давления и деформаций в виде интегральных выражений. Полученные 
решения позволяют рассматривать задачу с разными типами смешанных граничных условий: только по упру- 
ГОСТИ, ТОЛЬКО ПО фильтрации или с изменением типа обоих условий. 

Ключевые слова: пористость, осесимметричная задача, консолидация, неоднородное полупространство. 


Введение. Решению различных задач теории консолидации посвящено большое число работ. 
Краткий обзор основных работ опубликованных до 2001 года и посвящённых контактным задачам 
теории консолидации представлен в [1]. Большинство работ, в которых рассматриваются кон- 
тактные задачи для пороупругой среды, посвящены исследованию немонотонному изменению 
упругих свойств в покрытии. Целью данной работы является постановка осесимметричной задачи 
и представление её фундаментального решения для пороупругого полупространства с функцио- 
нально-градиентным покрытием. 
1. Основные уравнения теории деформаций пористо-упругих сред. 

Представим основные уравнения теории консолидации для пороупругого полупростран- 
ства. С полупространством свяжем цилиндрическую систему координат (г,ф,2). Обозначим че- 


рез и, и, ии смещения вдоль осей г,ф,2; 0,,0,,0,,Т„т — радиальное, угловое, нормаль- 


гф! Тр 
ное и тангенциальные напряжения соответственно. 
Для определения полей перемещений, деформаций, напряжений и порового давления в 
рамках квазистатической несвязанной задачи консолидации мы имеем следующие уравнения [1]: 
1. Определяющие соотношения связи напряжений и деформаций с поровым давлением 
возьмём в виде: 


_ 2Сбу 











[е] =+20Е. -ар, 
р И. 
2бу 
= ЕО, (1) 
В 2 =+2СЕ, -ар, 
1-2 
Т, =26Е„, Т, =2СЕ.,, Ть =26Е, 


1 Работа выполнена в рамках инициативной НИР. 
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3(% - у) 

В(1-2%)(1+%.) 

_ 268? (1 -2у)(1+у,} ‚ 268 (1+) 

9(%-у)(1-2м,) 3(1-2%) 

нированный коэффициенты Пуассона у<у, <0,5, Е=Е, + ВоЧЕЕЫ = деформация скелета, 
с — изменение содержания жидкости в единице объёма, В — коэффициент Скемптона. 
выд ВВ 
В-В. +т (В, -В.) 


где а= 





коэффициент эффективного напряжения Био, 





Е, у, У — соответственно дренированный и недре- 





‚ В. — сжимаемость твёрдой фазы, скелета; В, — сжимаемость поро- 


3(1 -2%) 


вого флюида, 1/7 — пористость, ах 
ыы . 6+) 


— дренированная сжимаемость скелета [3], 


С — модуль сдвига. 
2. Уравнения равновесия неоднородного по глубине полупространства при отсутствии 
массовых сил, записанные в цилиндрической системе координат (г, ф, 2: 
00, 1 0Тю , ОТ +. 9-95 _ 
ог Г 9 07 Г 
от 106, ОТ 2т 
ГР + 0, (2) 
ог Г 0 02 Га 
ОТ» + 10. + 00, + Те 
ог Г 0 02 Га 
Уравнения (1) и (2) определяют напряжённо-деформируемое состояние скелета и распределение 
порового давления. 
3. Движение поровой жидкости предполагается подчиняющимся закону Дарси. Закон Дар- 
си для квазистатического течения жидкости возьмём в виде: 


А о 
где к — коэффициент гидравлической проводимости (фильтрации), к=К /и>0 — коэффициент 
проницаемости, и — динамическая вязкость жидкости. 
4. Уравнение неразрывности квазистатического потока жидкости имеет вид 
06 оф 10, ф ду, 
+ 
Для определения положительной энергии деформации будем считать, что параметры имеют сле- 
дующие ограничения: С >0, 0<В <1, -1<у<\, <0,5, к>0. 
Наряду с модулем сдвига с и коэффициентом Пуассона у, для описания упругого пове- 


дения твёрдого изотропного тела используются коэффициенты Ламе Л и М или модуль Юнга Е и 
коэффициент Пуассона у. При приложении нагрузки, мгновенное деформированное состояние 
среды даётся решением задачи теории упругости с постоянными Ламе М и Л, =Л+а*и, то есть с 
модулем сдвига Си коэффициентом Пуассона 
Л-+а’н 
У Е р) ы (5) 
2 (^ +М-+а р) 

Коэффициенты М (иногда обозначаемый С и называемый модулем сдвига) и Ламе Л связаны с 


модулем Юнга Е и коэффициентом Пуассона у соотношениями 
И К Еу | МИА) ее (6) 
2(1+у) (1+у)(1-2у) М+Л 2(М+л) 





0, 

















= 0. (4) 
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Компоненты малой деформации выражаются через смещения по формулам Коши следующим об- 
разом: 





ди 1[ди ди 
ви +И|, Е, =—, 
г ? г| д 07 


1[ ди ду и 
= + | 
т? 2\гоф д г 


з кн" а кои. 
"262 0) 2102 Раф’ 
ди ду и ди 
=—+ +—+ : 
ОГ ГОф Г 02 
2. Постановка контактной задачи для градиентного пороупругого полупространства. 
Рассмотрим пороупругое полупространство, упругие характеристики которого непрерывно 
меняются с глубиной в пределах прилегающего к поверхности слоя толщины Н/ , а затем стабили- 
зируются и остаются постоянными. Пусть 2а — диаметр круга в основании жёсткого штампа. 
Предположим, что внутри круга приложена произвольная равномерная нормальная нагрузка (г) 
и касательные нагрузки &(7) и Е). Остальная часть поверхности свободна от напряжений. 
При осевой симметрии напряжённого состояния смещения, деформации и напряжения не 
зависят от угловой координаты ф. Из системы уравнений (2) в этом случае имеем 


09, ОТ 0,=0. 


Г [2 





(7) 

















ОР р г -° 

1. 

ОТ 00, т» (8) 

+ + = 0. 
ОГ 02 Г 
от от 2т 
2. + -=0. (9) 

[6] а 02 Г 


Система (8) описывает осесимметричное напряжённое состояние, возникающее, под действием 
нормальной к поверхности нагрузки, а уравнение (9) — равновесие полупространства, скручива- 
емого касательным усилием. 

Выражения для компонент деформаций (7) при этом упрощаются 








ди и ди/ 
= ' и = ' 
ог г 02 
1[0ди ди 1 ди 1[диу и 
=. = + ‚Е = , Е 5 й 10 
ы 25 я ТО, 55 :] то 
ди и ди 
= ++ —. 
ОГ Г 02 


Представления для напряжений (1) через смещения, с учётом (10) и (5) примут вид 
ди и ди 
в, =[^Л(2)+2М(2)}-—+Л(2)-+—|-а(2)р, 
‚Аааа име ар 


0 (^(2) +2М(2)) + ^(2) 
(11) 


в, = (А(2)+2м(=)) +2 


) 
а м2 + от =М(2 


02 0Г 











ду а(2)= 3(%-у(2)) 


То -М(2)- ы 
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В осесимметричном случае уравнение квазистатического потока жидкости (4) примет вид: 
И Е, (12) 
ог Г 02 
или, иначе, подставляя в (11) значения из (3): 
ь 2 к (2 
Ра (2) р _ 0, (13) 
0г ГОГ 02 02 02 








Подставляя выражения (11) в уравнения равновесия (8), получаем систему дифференциальных 
уравнений в частных производных второго порядка относительно перемещений и порового дав- 
ления. Эти уравнения могут быть представлены в форме 
и 0Е ди ди ор 
М(2)| Уи -— |+ (М(2)+Л(2))—+мМ'(2 + =а(2)—, 
(еды - + (мела) [5+ (2) 
М(2) Уи + (м(2) + ^(=)) = + 2м (=) + ^№'(2)==а'(2)р+а(2)р’, 























(14) 
у р ду 
м2 [9 - м (#5 о, 
к(2)У’р+к’(2)р’=0. 
Здесь используются следующие обозначения: 
М Л 2 
м(2)=2 (=) ^(2)=° и а ИЕ ; 
[674 Гор гог\ дг) 02? (15) 
сии, 0 ыы ан) р’_бР 
Ве @ ' а ' 02 


Мы полагаем, что через участок поверхности (круга радиуса а) происходит инфильтрация (отка- 
чивание) жидкости, остальная поверхность абсолютно непроницаема. Тогда, на поверхности гра- 
ничные условия имеют вид 


2=0: р’(г,0) =-В(г)к, г<а, 
р’(г,0)=0, гра, 
о, (г,0)=-9 (г), О<г<а, 
о, (г,0)=0, а<г<о, 
т» (г,0)=-а (г), О<г<а, те 
т» (г,0)=0, а<г<о, 
тб =Е (А). бегаа; 
То (^,0) =0, а<г<о 


На границе сцепления неоднородного слоя с однородным полупространством, при 2 =-Н , в силу 
непрерывности, должны выполняться условия сопряжения по смещениям, напряжениям и поро- 
вому давлению. 


об (г,-Н)= с° (г,-Н), < (г,-Н)= т” (г,-Н), 
иб (г,-Н) = и? (г,-Н) им (г,-Н)=и? (г,-Н), (17) 
рб (",‚-Н)= р (",‚-Н), (2°) (",-Н)= (2°) (",-Н). 
На бесконечности, при (г,-2)-› смещения, деформации и напряжения исчезают. Значение 


порового давления при этом постоянно. 
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ИИ (МЕ, ЕЕ, 8„,0,,быб„,т„,р’)=0, 


"пр-р., Итр=р. (18) 


2——0 
Ия (МИ, Е.Е, Е,,Е Об бытыР О: 


171 


Учитывая, что в пределах прилегающего к поверхности слоя толщины Н характеристики рас- 
сматриваемого полупространства непрерывно меняются с глубиной, полагаем, что коэффициент 
Пуассона у(2), модуль сдвига С (2), коэффициента Скемптона В(2) и коэффициент фильтра- 


ции к(2) являются непрерывными функциями координаты 2 , такими, что 


1. 6(2)=6° (2), у(2)=\ (2), 
к(2)=к< (2), В(2)=В° (2), -Н<2<0 
2.6(2)=6(-Н)=6°, у(2)=у(-Н) =, (19) 
к(2)=к(-Н)=к®, В (2) =В (-Н) = В®, —ю<2<-Н 
3. 6° =6° (-Н), у =\ -Н), 


к =кб (-Н), В® =В< (-Н). 
Индекс $ соответствует подстилающему однородному полупространству, а С — неоднородному 
слою. 
3. Построение фундаментального решения квазистатической осесимметричной задачи 
консолидации для неоднородного по глубине полупространства. 
Будем разыскивать решение для смещений и,и/ и порового давления р в виде интегра- 
лов Ханкеля 


) 


[и ‚(уг)уау, и(г,2) М (\у,2) 7, (уг)уау, 
о (20) 


РИ (уг)уау, р(г,2) ИЕ (у,2) 7. (уг) уау. 


Подставив (20) в систему дифференциальных уравнений в частных производных (14) и, прирав- 
няв к нулю подынтегральные выражения, получим систему обыкновенных дифференциальных 
уравнений второго порядка (здесь ' указывает на дифференцирование по 2 ). 
МИ" +у(М+л)и/" - у? (2м+л)И +МИ' + УМИ/ =ауР, 
(2м+^)и/"-у(2мМ+ Л)” - У?МИ/ + (2мМ' + Л’) И/' -УЛ'И =а'Р чар’, 


(21) 
МИ" +М'И ' - у?МИ =0, 
к.Р”-уж.Р+к’.Р’=0 
Граничные условия (16) примут следующий вид: 
(2м(0) +л(0)}и/'(у,0) - ул (6) (у,0)- а(6)Р (\,0) =-9 (у), 
м(о)(уи (у,0) + (\+0)) = В (у), м (6)и'(у,0) =Т (У), 
Р’(у,0) =-к-У (у), 
а а (22) 
@(у) = [9 (2). (ру)рар, Б(у)= [а (р).7 (ру)рар, 
(у) = [В(2)7 2. (ру)рар, Т (у о ‚ (ру)рар. 


Третье уравнение системы (21) не зависит от составляющих остальных уравнений этой системы. 
Его решение представлено в [4]. 
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Система уравнений, описывающая напряжённое состояние пористой среды под воздействием 
нормальной и радиальной нагрузках имеет вид 


МИ" +у(М+л)и/' у? (2м+ Л) И +М" + УМ'И/ = ауР, 
(2м+л)иИ/"-у(2м+ ^)И' - У?МИ/ + (2М' + Л’) И/' - УЛ’ =а'Р+аР', (23) 
к(2).Р"-уж(2)Р + (к) (2)-Р'=0. 
С граничными условиями 
(2м(0)+л(0)} и/'(у,0)-у ни __ в(0)Р (\,0) = -@ (у), 
ом (о) у, — р (у), 


Р'(у, _ К 7 (24) 


м Л (ру)рар, В(у )= [92 ‚ (ру)рар, 


)= В (р) 7» (ру)рар. 


Решение аналогичной системы для случая о деформации представлены в [5]. По- 

скольку уравнения теории термоупругости по форме совпадают с уравнениями теории изотерми- 

ческой деформации для пористой среды [6], содержащей вязкую сжимаемую жидкость, то даль- 

нейшие рассуждения проведём в соответствии с методом решения, предложенным в работе [5]. 
Введём обозначения: 


х ах х,,х ), 
х, =Ц, х, =И' х.=И/, Х.=И/', (25) 
Х; =Р, Хх, =Р" 


Запишем систему (23) в матричном виде, при этом явно выделим части соответствующие покры- 
тию и подложке: 









































С 
дсхе, -Н<2<0. (26) 
[674 
0 1 0 0 0 
52М+Л м м М+л ца 
‘м М Ум м М ) 
0 0 0 1 0 0 
АС = ^ М+л о), м 2м + а’ а |; 
мал У2м+л "2М+л 2М+Л 2М+Л ЭМ+Л 
0 0 0 0 0 1 
0 0 0 0 у? се 
К 
5 
К АЗх5, фо 2<-Н. (27) 
[674 
0 1 0 0 0 0 
22М+Л 0 _,М+А уа 0 
М м м 
д5_| 0 0 0 1 оо | 
о М+Л у? М о о а 
2М+ Л 2М+ Л 2М+ Л 
0 0 0 0 0 1 
0 0 0 0 у? 0 
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Граничные условия при этом имеют вид 
(2м(0) +^(0)) хх (у,0) - Ул (0) хх (у,0) -а(0)х; =-© (у), 
м(0) (уж (%0) + ) 
) У), (28) 
с (%,0) = х? (у,-Н), ху (у,0) х 
(у,0) ыы № (у,-Н), № (у,0) — 4 (у,-Н), 
х5 (\,0) =хё (у,-Н), хе (у,0) = хе (у,-Н). 
Общее решение системы (23) для однородного полупространства Л’ =М' =а'=к’ = 0, 
М> 0, Л>0, а>0, к>0 имеет вид 


Хх 
Хх 


х? (\,2) = |4 + 24, . |= (к, т ве", 
х (У,2) -[^ + (1 * у2) 4. + =. (к. + ке ме", 
хз (\,2) = |4 + (-к, +\2) 4, +5 (к, + ке", (29) 


хе 6,2) [а + (1-к, +2) а, + [+ на] ме", 


= вах уе, 
_ Л+3М К К 


СРзм' 74 мМ+л) 37 4м 


где а, (/=1,2,3) — произвольная функция параметра у. 





Решение хб (\,2) системы дифференциальных уравнений (26) строится методом модули- 
рующих функций [7]. Будем искать х° (у,2) в виде 
х (,,2) => (У)а, (\,2)е*. (30) 
Векторы а, (\,2), (/=1,2,3) определяются из решения следующей задачи Коши: 
ей = Аба -уа, -Н<2<0, /=1,2,3. (31) 


[74 
при начальных условиях для 2 =-Н 


а, азы = (1,1, у,0,0) 
а, И = (ух, у+ 22, к, +у2,у-ку +72, 0,0) | 





2=- 





а, (\,2)|,_н = [= (ку +к,)2,к. + (к, + коде, 2+ (к, +к,)2,к, + (к. + омеди 
2=-Н 


Константы а, (у) (/=1,2,3) определяются из условия (32). Таким образом, мы имеем 


ха, (ум, ()=В(у), (32) 
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м, (у) =-^(0) уа; (у,0) + (^ а. )) а’ (\,0)-а(0) а? (у,0) 
Е. ) =М(0)а, (у, о р уа; (у,0), 
0, (у) = а° (ч.0), 0, (у)= 0, (у)=0. 
где а’ (\,2), (!=1,2,3) обозначает А-ю компоненту вектора а, (у,2), / =1,2,3, К =1,2,3,4,5,6. 
а, =- К, (у)® (у) № (у) (у) К» (у)к-У (У), 
д, - ко) р (чеку), (33) 
Во 


























: А,» (У) . Л, (\)" " д, (У) О: (У) 
ка Ку Ка 


А, (У) = М, (У)М, (у) - М, (У)М, (у). 
Окончательно получаем следующее выражение для компонент вектора решения хс (\у,2) при 
2>-Н 





И(у,2) = х, = (4 (у,2)9 (у) В (у,2)Б (у) + В (у,2)ку (у))е* м, 
И" (№2) =х, = (В (у,2)9 (у) +6 (у,2)В (у) +В (у,2)ку (у))е" №, 
ИИ (\,2) = х, = (В (\,2)9 (У) +6 (у,2)В (у) + В (у,2)кУ (у))е" м, 
И/'(у,2) = х+ = (В (у,2)0 (у) +В (у,2)Б (у) + В (\,2)ку (у))е* № 
Р(у,2)=х, = В (у,2)кУ (у)е*/у, (34) 
Р'(у,2) = Хз = [8 (у,2)кУ (у)е\ /\ 
р (у,2) = (К (\)а; (у,2) - ® (уда, (у,2))\, 
ь (2) = (К (у)а; (№2) К (у)аи (у,2 у, 
ибо даи)- воаци- 


Введём обозначения 


1 (г.г) = [1 (ч,2)9 (уе (уг)ау (к = ь2 11) (К =3,4 1 =0) 








= [№ (у,2) (ед (мг) (К =1,2 11) (К =3,41=0) 











14 (", ик 2) (\)е* 2 (уг)ау, (к=1,21=1),(К=3,4,5,61=0) 





(35) 








[и (у)е*“ 4, (уг)уау, (К =21=0),(К =47=1) 








р(уе* 1 (уг) уу), (к=21=0), (К =41=1) 








Ще (уУ)е*, (уг)уау, (К =21=0),(К =41=1) 
0 
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В соответствии с (10)—12) запишем выражения для смещений, порового давления и деформаций. 
3 3 
и(г,2) =-У 1, (г,2) им (г,2)= УТ (г,2), 
1=1 1=1 


Ор (г,2 
ВЕ), Е 
3 


} 1 (Ре) 








ы (36) 








№ 
ей 


Выражения для напряжении принимают вид 
с, (ме) м)- Ди ("2+ 
=1 [ 
ме -ЕУ (+) (г,2)|-а(2)Ть (г,2), 


(37) 


= 
ты 
2 
5 
- 
х 
м 





о т а 


Заключение. В работе рассмотрена постановка осесимметричной задачи для пороупругого по- 
лупространства с функционально-градиентным покрытием. Получены представления для напря- 
жений, смещений, порового давления и деформаций, возникающих под воздействием равномер- 
ной нормальной и касательных нагрузках. 
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АХТ$УММЕТВТС СОМТАСТ СОМ$ОЕТОАТТОМ РВОВЕЕМ РОВ СОМТТМООЦ$УЕУ 
МОМНОМОСЕМЕОЧЦ$ ТМ ОЕРТН НАЕЕ-5РАСЕ* 


Е. М. Уемсь 
(Роп Зае Тесптка! Упмег$Ку) 


Тре чейетепЕ о! те ахвутте тс сопзо!аайоп рго Мет Юг а рогоиз, поп-потодепеоиз т-дерЁй Ва! -5расе, апа 
те соп$гисйоп оЁ пе паатета!/ оивоп ю аейптд {те @расетепЕ Пе/45, #а/п5, ге55ез, апа роге ргез- 
5иге ипаег пе дмеп Боипаагу сопаопз аге сопз!аЕгеа, ТПе 5о/иНоп ю {Те едиаНоп$ аезс"Ытд {те рогоиз те- 
ит $#те5 ипдЕг {те реаг Гогсе, апа ипдЕг те погта! апа габ!а! [сасб 15 сопу!егеа абзо/ше/у апа йгезресвуе- 
у Напке! пиедга! {гапюЮгт иНсВ аШои/5 сопуегта {те ргоШет {0 {те зо/иНоп о! {пе огтагу тегепва! зесопа- 
ОгаЕг зузет 15 изеа. Ехргезюп$ Юг $#ез5ез, @ср/асететЁ, роге ргезиге, апа тат т {те ютт оЁ п®едга! 
ехрге5зюп5 аге обеатеа {гоидй пе ятиайпа Гипсйоп те од. Тре гезийта зоийоп5 аШои/ сопуаетта пе 
роет ийн уапоиз {урез оЁ {те тихеа Боипаагу сопа@!Нопз: оту оп еазНсйу, оту оп ВЁгайоп, ог ий Ве 
спапде оЁ те ёуре о! Бо сопа!опз. 

Кеуигогаб&: рогозйу, ахутте "с ргоМет, сопзо!аавоп, поп-ротодепеоих Ве! 5зрасе. 


1 Тпе гезеагсн {5 допе мт {пе Нате о! {Не паерепдепе В&0. 
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